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HIGHLIGHTS 


►  Conduct  a  unique  in  situ  AFM  measurement  on  micro-machined  sample. 

►  Exclude  secondary/non-active  materials  and  remove  vagueness  due  to  irregular  geometry. 

►  Develop  and  validate  a  proposed  coupled  electrochemical  and  mechanical  model. 

►  Reveal  the  key  effects  of  geometry  and  kinetics  on  active  material  stability. 
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Structural  instability  due  to  intercalation-induced  stresses  in  electrode  materials  is  one  of  the  key 
degradation  mechanisms  of  Li-ion  batteries.  Fragmentation  of  material  degrades  structural  integrity  and 
electrical  resistance,  and  also  accelerates  harmful  side  reactions.  In  situ  experiments  are  the  appropriate 
approach  for  investigating  the  actual  time-dependent  nature  of  the  behavior  changes  of  an  electrode 
material  while  it  is  charged  and  discharged.  In  the  current  work,  a  unique  in  situ  electrochemical  atomic 
force  microscopy  (ECAFM)  measurement  is  made  on  samples  of  cylindrical  shape,  which  are  micro- 
machined  by  focused  ion  beam  (FIB)  microscopy.  This  pre-defined  geometry  allows  the  exclusion  of 
secondary,  non-active  materials  from  the  electrochemically  active  material  as  well  as  the  removal  of  any 
vagueness  owing  from  the  irregular  geometry  of  particles.  The  experimental  results  are  also  used  to 
validate  a  proposed  coupled  electrochemical  and  mechanical  model  for  determining  the  stress— strain 
state  of  active  electrode  material  during  electrochemical  cycling.  The  results  produced  using  this 
model  correlate  strongly  with  the  experimental  data.  The  combined  results  reveal  the  key  effects  of  the 
geometry,  kinetics,  and  mechanics  of  electrode  materials  on  the  stress— strain  state,  which  acts  as 
a  barometer  of  the  structural  stability  of  a  material. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

L1C0O2  in  bulk  powder  form  has  been  commonly  used  as 
a  cathode  material  for  commercial  Li-ion  batteries  in  consumer 
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electronics  since  1990  [1].  For  increased  safety  in  large-scale 
batteries  for  electric  vehicles  and  grid-scale  storage,  other  deri- 
vates  of  LiCo02  are  used  [2,3].  Due  to  this  history  in  consumer 
electronics  and  the  large  body  of  associated  research,  the  electro¬ 
chemical  behavior  of  LiCoCb  is  well  known  and  therefore  a  perfect 
candidate  for  developing  and  testing  new  characterization  methods 
in  association  with  theory.  We  selected  this  material  for  charac¬ 
terization  by  a  new  electrochemical  atomic  force  microscopy 
method. 

Thin  films  for  use  as  cathode  materials  have  been  fabricated  by 
diverse  methods  including  e-beam  evaporation  [4],  pulsed  laser 
ablation  [5],  sputtering  [6],  and  electrostatic  spray  deposition  [7]. 
Generally,  upon  deposition  of  the  thin  film,  the  cathode  is  annealed 
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above  700  °C  to  improve  its  crystallinity  [8-10].  The  UCo02  thin 
film  has  a  rhombohedral  structure  consisting  of  layers  of  close- 
packed  networks  of  oxygen  ions  separated  by  alternating  layers 
of  Li  and  Co  ions.  The  hexagonal  cell  has  lattice  parameters  of 
a  =  2.815  A  and  c  =  14.05  A  [11]. 

During  cell  cycling,  anisotropic  structural  changes  occur  as  Li- 
ions  (de intercalate  through  the  active  materials.  X-ray  diffraction 
measurement  has  observed  the  change  of  lattice  spacing  in  LixCo02 
as  a  function  of  x,  where  the  change  in  the  c-axis  is  much  greater 
than  that  in  the  a-axis.  The  spacing  change  in  the  c-axis  increases 
during  charging,  reaching  a  maximum  at  around  x  =  0.5,  after 
which  it  decreases.  Increasing  electrostatic  repulsion  between  Co02 
layers  due  to  deintercalation  causes  the  c-axis  lattice  parameter  to 
expand  by  2.4%  from  c  =  14.08  A  in  Lio.9gCoi.oi02  to  c  =  14.42  A  in 
Lio.4gCoi.oi02;  the  a-axis  only  shrinks  by  0.34%  from  a  =  2.816  A  in 
Li0.ggCoi.oi02  to  a  =  2.807  A  in  Li0.4gCoi.oi02  [12].  These  changes 
have  been  confirmed  by  several  in  situ  and  ex  situ  X-ray  powder 
diffraction  experiments  [13-15]. 

Structural  instability  in  the  form  of  micro-cracks  caused  by 
cycling-induced  changes  in  the  lattice  parameters  has  been 
proposed  as  a  major  cause  of  performance  degradation  in  LiCo02 
batteries.  These  micro-cracks  can  lead  to  fragmentation  of  the 
particles,  which  can  affect  the  particle  network  and  electrical 
resistance,  as  well  as  accelerate  side  reactions  with  the  electrolyte. 
A  TEM  examination  of  LiCo02  cathodes  cycled  between  2.5  V  and 
4.35  against  a  Li  anode  [16]  has  shown  that  a  significant  fraction  of 
the  particles  was  severely  strained,  fractured,  and  possessive  of 
a  high  density  of  extended  defects.  Similar  damage  may  be 
responsible  for  degradation  in  the  long-term  cycling  of  LiCo02  Li- 
ion  batteries.  Dislocations  in  LiCo02  were  observed  by  trans¬ 
mission  electron  microscopy  [17].  It  was  found  that  the  configu¬ 
rations  of  all  dislocations  are  glissile  and  indicative  of  active  slip 
planes.  These  glissile  dislocations  might  suggest  possible  damage 
mechanisms  during  cycling. 

While  pulverized  particles  have  been  observed  in  cycled  mate¬ 
rial,  investigation  into  changes  in  stress-strain  state  and  damage 
accumulation  as  a  time-dependent  phenomenon  is  still  under  way. 
Experimentally,  such  works  require  in  situ  techniques,  which 
impose  certain  challenges  on  the  design  of  the  experimental  setup. 
The  Acoustic  Emission  (AE)  method  has  been  employed  to  register 
acoustic  waves  emitted  from  the  composite  Si  electrode  when 
fracture  occurs  during  charge  and  discharge  [18].  Recently,  nano¬ 
mechanical  quantification  of  the  elastic,  plastic,  and  fracture 
properties  of  LiCo02  was  examined  by  measuring  the  Young’s 
elastic  modulus,  hardness,  and  fracture  toughness  [19].  Further, 
design  criteria  for  electrochemical,  shock-resistant  battery  elec¬ 
trode  were  made  based  on  the  analytical  model  and  in  situ  acoustic 
emission  experiments  [20]. 

While  this  technique  allows  for  the  quantification  of  freshly 
produced  cracks  at  each  particular  time  interval,  the  measurement 
of  strains  in  the  material  with  the  purpose  of  understanding  the 
constitutive  electrochemical— mechanical  behavior  is  limited.  In 
situ  X-ray  experiments  provide  information  on  strain  and  phase 
transformations  in  cathode  materials  [21-23];  however,  the  tech¬ 
nique  focuses  on  the  atomistic  scale  and  cannot  give  a  detailed 
analysis  of  the  electrode  scale.  Use  of  Scanning  Probe  Microscopy 
(SPM)  techniques,  on  the  other  hand,  enables  the  observation  and 
collection  of  data  with  rather  high  resolution  on  a  nano-meter 
scale.  While  early  SPM  experiments  targeted  observation  of  the 
overall  roughness  change  in  the  electrode  during  electrochemical 
cycling  [24],  the  technique  has  been  further  developed  to  probe 
diffusion-induced  strains  in  small  features  on  thin  film  Li  ion 
battery  cathodes  [25,26].  In  situ  AFM  measurement  for  roughness  of 
sputter-deposited  Si-Sn  thin  films  for  electrode  material  has  been 
conducted  [24].  An  in  situ  electrochemical  AFM  measurement  was 


conducted  to  study  surface  and  dimensional  changes  of  individual 
LixCo02  crystals  during  lithium  deintercalation  [27].  The  extent  of 
expansion  measured  along  the  Chex  axis  corresponded  with 
previous  X-ray  powder  diffraction  results. 

Though  previous  work  has  attempted  to  monitor  changes  in 
lattice  parameters  or  surface  roughness  during  intercalation/dein¬ 
tercalation,  the  results  were  obscured  by  the  effects  of  particle 
geometry  and  particle  size.  As  previously  mentioned,  lattice 
expansion  in  LixCo02  is  a  function  of  state  of  charge  [28],  which 
means  that  stresses  due  to  volume  changes  result  from  a  concen¬ 
tration  of  Li-ions  in  the  cathode  particles.  Depending  on  the 
geometry  and  the  size  of  the  particles  being  measured,  the 
concentration  profile  changes  significantly.  Furthermore,  a  bulk 
composite  cathode  includes  additive  particles  such  as  carbon  black 
particles  and  PVdF  binders  [28].  These  additive  materials  may  affect 
the  stress  level  in  active  particles  by  changing  the  boundary 
conditions.  To  date,  there  has  been  no  attempt  to  measure  the 
strain  on  an  electrode  of  pre-defined  geometry  where  the  shape 
and  volume  changes  can  be  observed  and  measured  accurately 
without  being  hindered  by  the  presence  of  additive  materials. 

In  addition,  it  is  nearly  impossible  to  maintain  a  sample  in 
a  perfectly  level  position  during  scanned  probe  measurements  on 
nanometer-scale  geometries.  Even  a  small  amount  of  tilt  can 
adversely  affect  the  accuracy  of  the  results.  Also,  the  AFM  image  is 
a  result  of  the  interaction  of  the  tip  shape  with  the  surface  topog¬ 
raphy;  the  curvature  and  tip  sidewall  angles  affect  the  image.  For 
these  reasons,  the  obtained  image  may  need  to  be  processed  for  the 
sake  of  better  visual  effect  or  to  eliminate  measurement  error.  Here, 
a  well-defined  geometry  is  a  benefit  because  correction  is  easier 
than  in  the  case  of  randomly-shaped  geometry. 

Modeling  efforts  on  a  particle  level  have  been  based  on  the 
analogy  between  thermal  stresses  and  stresses  introduced  by  the 
diffusion  of  a  chemical  species  into  a  material  [29-35].  A  model 
utilizing  binary  diffusion  coefficients  has  been  developed  and  used  to 
compute  the  stresses  in  LiMn204  [36]  and  carbonaceous  material 
[37].  The  thermal  analogy-based  modeling  has  been  extended  from 
a  single  particle  to  a  composite  electrode  [38],  and  the  stress  state 
produced  by  lithiation  was  used  as  an  input  for  the  fracture 
mechanics  analysis  of  an  existing  crack  in  active  material  particles 
[39,40].  The  description  of  volumetric  expansion  due  to  intercalation 
based  on  equations  of  thermal  expansion  is  rather  appealing  due  to 
the  simplicity  of  the  elastic  equations,  and  this  approach  is  adopted  in 
the  current  investigation  by  considering  electrochemical  interactions 
between  the  electrode  micro-specimen  and  the  electrolyte. 

In  the  present  investigation,  a  micro-specimen  of  pre-defined 
geometry  was  fabricated  and  an  in  situ  AFM  technique  was 
applied  to  measure  the  morphology  changes  in  this  sample  during 
charging  and  discharging.  A  Focused  Ion  Beam  (FIB)  microscope 
was  used  to  mill  a  cylindrical  sample  out  of  a  LiCo02  thin  film 
electrode.  Simultaneously,  a  mechanical-electrochemical  model 
including  thermodynamics,  electrode  kinetics,  and  transport 
phenomena  was  developed;  this  model  was  used  to  simulate  the 
same  configuration  as  that  in  the  experiment. 

The  present  work  thus  attempts  to  confirm  the  mechanism  that 
correlates  structural  instability  with  capacity  fade  in  LixCo02  elec¬ 
trodes.  This  work  focuses  on  examining  the  surface  morphology 
evolution  of  a  LixCo02  thin  film  electrode  via  in  situ  AFM 
measurement  of  a  micro-electrode  and  microscale  electrochemical 
modeling.  In  this  study,  we  have  the  following  objectives: 

1.  Conduct  an  in  situ  AFM  study  to  measure  morphology  changes 
of  a  pre-defined  thin  film  based  micro-electrode  as  a  function 
of  state  of  charge,  specimen  size,  and  potential  sweep  rate. 

2.  Develop  a  multiphysics-based  model  for  calculating  the 
volume  expansion  and  deformation  due  to  intercalation 
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considering  the  multiphase  interaction  between  the  electrode 
and  the  electrolyte. 

3.  Compare  the  simulation  results  with  measurement  results  and 
further  extend  the  prediction  of  mechanical  characteristics  of 
the  specimen  during  (de  intercalation,  which  cannot  be 
measured. 

2.  Experiments 

2.1.  Synthesis  ofLiCo02  thin  film  electrode  and  micro-machining 
process  using  FIB 

Thin  film  samples  were  prepared  by  radio  frequency  magnetron 
sputtering  of  a  cold-pressed  and  sintered  UCo02  target.  The 
substrate  consisted  of  polished  99.6%  Alumina  with  a  300  nm  thick 
gold  film  deposited  as  a  current  collector.  In  order  to  improve 
adhesion  of  Au  to  the  substrate,  a  20  nm  thin  layer  of  Ni  was 
deposited  first.  UC0O2  was  deposited  at  an  average  deposition  rate 
of  18  A  min-1,  to  a  final  thickness  of  2  pm.  The  samples  were 
annealed  at  700  °C  for  2  h  under  a  flow  of  oxygen. 

Fabrication  of  the  micro-specimens  was  performed  in  a  dual 
beam-focused  ion  beam/electron  beam  microscope  (Hitachi 
NB5000).  The  described  procedure  yielded  thin  films  with  rather 
low  surface  roughness  and  small  grain  size  (about  50-70  nm)  of 
UC0O2.  As  a  result,  no  surface  deposition  of  carbon  or  tungsten  was 
necessary  in  order  to  achieve  smooth  surfaces  on  the  milled  pillars. 
In  order  to  minimize  the  implantation  of  the  Ga  ions  into  the 
surface  layer  of  the  material,  fabrication  was  performed  with 
decreasing  current  in  each  run  (ranging  from  3.79  nA  to  0.07  nA  in 
the  final  polishing  step).  The  fabrication  was  monitored  with  SEM 
in  order  to  make  sure  the  sample  remained  attached  to  the  current 
collector.  A  thin  film  with  a  FIB-fabricated  micro-pillar  produced  in 
this  manner  is  shown  in  Fig.  1. 

2.2.  Electrochemical  measurement 

The  voltage  profile  of  U/UC0O2  to  be  used  as  a  reference  in  in 
situ  AFM  measurements  was  obtained  from  a  thin  film  sample 
cycled  in  a  Swagelok  cell.  The  cell  was  assembled  in  an  argon-filled 
glove  box  with  H20  maintained  at  less  than  1  ppm.  The  cell  con¬ 
sisted  of  a  UCo02  thin  film  cathode  electrode,  a  lithium  foil  anode 
electrode,  and  a  polypropylene  separator  (Celgard  2500).  The 
electrolyte  solution  consisted  of  1  M  LiPF6  in  an  1:1  v/v  mixture  of 
ethylene  carbonate  (EC)  and  dimethyl  carbonate  (DMC).  The  Swa¬ 
gelok  cell  was  allowed  to  rest  overnight  in  the  glove  box  after 
assembly  prior  to  cycling.  The  initial  formation  cycle  was  conducted 


Fig.  1.  SEM  image  of  the  specimen  including  four  pillars.  A  larger  pillar  (no.  3, 
r0  =  3.8  pm)  and  a  smaller  pillar  (no.  2,  r0  =  2.05  pm)  are  selected  for  quick  scanning, 
and  are  marked  with  the  dotted  area.  One  AFM  example  is  shown  in  the  top-left 
corner. 


at  a  rate  of  C/10  between  3.0  and  4.5  V  versus  lithium  using  a  Bio- 
Logic  VMP3  potentiostat. 

2.3.  In  situ  AFM  measurement 

In  situ  AFM  experiments  for  the  morphology  evolution  of  the 
samples  were  conducted  with  a  Bruker-Nano  Multimode  AFM  with 
a  Nanoscope  III  controller  coupled  with  a  BioLogic  VMP3  cycler.  The 
AFM,  sitting  on  a  vibration  isolator,  was  located  inside  an  argon- 
filled  glove  box,  with  only  the  controller  and  the  potentiostat 
located  outside  the  glove  box.  The  half  cell  was  assembled  using 
a  Bruker-Nano  MMTEC  AFM  Fluid  cell  inside  the  glove  box. 

A  schematic  diagram  of  the  AFM  Fluid  cell  cross-sectional 
configuration  is  shown  in  Fig.  2.  Lithium  foil,  used  for  the  counter 
electrode,  was  placed  on  top  of  a  copper  current  collector.  The  AFM 
Fluid  cell  was  placed  on  top  of  the  working  electrode,  and  the  rim 
was  sealed  with  a  silicone  O-ring.  The  sample  was  placed  in  the 
center  below  the  silicon  nitride  AFM  scanning  tip  (SNL,  Bruker- 
Nano  Probes).  Finally,  the  electrolyte  was  injected  from  the  right 
inlet  port  with  the  excess  exiting  from  the  left  inlet  as  shown  in 
Fig.  2. 

To  study  the  time-dependent  morphology  changes  in  the 
particle,  cyclic  voltammetry  was  employed.  A  linearly  changing 
ramp  voltage  was  applied  to  the  electrode.  The  rate  of  potential 
change  (sweep  rate)  was  set  at  two  values,  0.39  mV  s-1  and 
0.78  mV  s-1.  While  the  half-cell  was  potentiodynamically 
controlled,  in  situ  AFM  images  were  collected  in  contact  mode 
with  a  non-conductive  silicon  nitride  tip.  The  stiffness  of  the  tip 
was  0.58  N  m-1  and  the  radius  was  10  nm.  Fig.  1  shows  an  SEM 
image  of  a  specimen.  This  specimen  has  four  pillars  and  one 
reference  at  the  center  of  the  specimen  where  the  LiCo02  thin  film 
was  removed  by  FIB,  exposing  the  Alumina  substrate.  To  capture 
the  volume  change  of  the  sample  at  a  specific  voltage,  the  scan¬ 
ning  was  done  quickly  over  a  small  area.  It  was  found  that  side 
reactions  such  as  gas  generation  and  oxidative  decomposition  of 
the  electrolyte  significantly  reduce  the  image  quality,  further 
necessitating  quick  scanning.  Consequently,  the  minimum  area 
required  to  make  visible  the  effects  of  shape  change  was  selected. 
The  maximum  scan  height  was  5  pm  in  the  z-direction,  and  the 
scanning  area  was  8.75  pm  x  35  pm.  The  previously  mentioned 
reference  point  on  the  samples  was  included  to  correct  for  any 
drift  that  might  occur  in  the  sample  image  during  data  acquisition. 
Taking  all  limitations  into  account,  the  area  selected  for  AFM 
imaging  included  pillar  no.  3  (3.8  pm  in  radius),  pillar  no.  2 
(2.05  pm  in  radius),  and  the  reference,  as  marked  with  the  dotted 
rectangular  box  in  Fig.  1.  One  AFM  image  example  is  shown  at  the 
top-left  in  Fig.  1. 


AFM  tip 

electrolyte  inlet 


electrolyte  j  sample 
substrate 


Fig.  2.  A  schematic  diagram  of  the  AFM  fluid  cell  cross-section  configuration.  The 
LiCo02  thin  film  electrode  sits  on  a  bottom-insulated  substrate.  An  aluminum  current 
collector  is  located  between  the  thin  film  and  the  substrate.  Lithium  metal  wire  is 
attached  to  a  copper  current  collector. 
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2.4.  Image  post-processing 

The  obtained  image  was  processed  to  correct  for  the  tilt  of  the 
sample,  which  appeared  to  be  inevitable  and  might  have  intro¬ 
duced  errors  if  not  appropriately  treated.  As  can  be  seen  in  Fig.  3, 
the  tilted  plane  has  an  angle  of  6  with  the  horizontal  plane.  For  the 
correction,  we  assumed  that  the  plane  was  tilted  down  and  away 
from  the  upper  right  corner  of  the  specimen  as  viewed  in  Fig.  3.  The 
constants  in  the  equation  for  the  tilted  plane,  z  =  ax  +  by  +  c,  were 
found  based  on  least-mean-squares  fitting,  after  which  the  cross 
line  and  the  angle  between  two  planes  (z  =  ax  +  by  +  c  and  z  =  0) 
were  calculated.  In  order  to  perform  the  correction,  the  tilted  plane 
was  rotated  around  the  cross  line  between  the  two  planes.  After 
rotating  the  tilted  plane  through  the  cross  line  passing  a  point  A 
with  a  direction  vector  U  by  an  angle  6,  the  new  position,  X'  (located 
in  the  horizontal  plane  of  a  point  X  in  the  tilted  plane),  can  be 
calculated  from  the  product  of  sequential  rotation  matrices: 


*  all  lengths  are  in 


a  before  correction:  3D  view  C  after  correction:  3D  view 
r, - - - - - - -  c 


'0  10  20  30 

b  before:  side  view 


0  10  20  30 

d  after:  side  view 


X  = 


A  + 


(XA>uu 

|U|2 


X  (X-  A)smd 


(1  -  cos  6)  +  cos  0X  +  —U 


(1) 


Fig.  4  shows  one  example  of  a  post-processed  result  where  (a)  is 
the  obtained  image  and  (b)  is  the  cross  section  of  the  obtained 
image  showing  the  tilt.  Fig.  4c  and  d  shows  the  corrected  image 
after  rotating  the  image,  where  the  tilt  has  been  corrected. 


3.  Numerical  simulations 


Fig.  4.  One  example  of  post-processing.  The  original  image  and  corrected  images  are 
shown  as  indicated:  (a)  the  original  3D  image,  (b)  the  original  cross  section  image,  (c) 
the  corrected  3D  image,  and  (d)  the  corrected  cross  section  image. 

Similarly,  the  redox  phenomenon  occurs  at  the  interface 
between  LixCo02  and  the  electrolyte  as  expressed  below: 

discharge 

Li1_xCo02  +xLi+  -\-xe~  -<  ^  LiCo02  (3) 

charge 

The  rate  of  the  reactions  represented  by  Eqs.  (1)  and  (2)  can  be 
modeled  by  the  Butler— Volmer  equation: 


3.1.  Description  of  the  modeling 


Jt‘  =  lo[exP[^fvi )  -  exp 


(4) 


Mapping  of  the  Li-ion  concentration  inside  the  electrode  micro¬ 
sample  is  necessary  for  tracking  the  morphology  change  of  the 
sample  and  for  determining  strain  in  the  micro-pillar.  In  addition, 
the  Li-ion  concentration  in  the  electrolyte  around  the  micro-pillar 
is  needed  to  determine  the  kinetics  associated  with  the  reaction 
at  the  interface  between  the  specimen  and  the  electrolyte.  This 
means  that  a  diffusion-only  model  of  the  Li-ion  in  the  solid  particle, 
as  assumed  in  the  previous  work  [41  ],  is  insufficient  for  capturing 
the  exact  phenomena  at  the  interface  between  the  particle  and  the 
electrolyte.  At  the  interface  of  the  lithium  metal  and  the  electrolyte, 
Li-ions  are  dissociated  from  or  reduced  into  the  lithium  metal 
during  the  discharge-charge  processes  as: 


where  z'o  is  the  exchange  current  density,  aa,  ac  are  anodic  and 
cathodic  transfer  coefficients  for  the  electrode,  respectively,  R  is  the 
gas  constant,  F  is  the  Faraday’s  constant,  and  T  is  the  temperature. 
Here,  rp  stands  for  the  surface  overpotential  defined  as 

Vi  =  4>i  —  4>e  —  Uj  (5) 

where  (ph  (pe  are  the  surface  potentials  at  the  interface  of  the  electrode 
(i  =  1  for  the  anode  and  i  =  2  for  the  cathode)  and  the  electrolyte,  and 
Ui  is  the  open  circuit  potential.  The  exchange  current  density  deter¬ 
mines  the  rate  of  the  reaction  and  depends  upon  the  reactant  and 
product  concentration  adjacent  to  the  interface  as  follows: 


discharge 

Li  <  ^  Li+  +  e-  (2) 

charge 


io  =  k(csrc(c™x-Cs)a>(ce)a>  (6) 

where  cs  and  ce  represent  Li-ion  concentration  at  the  particle 
surface  and  in  the  electrolyte,  c™ax  stands  for  the  maximum  surface 
concentration,  and  k  is  the  reaction  rate  constant. 

The  governing  equations  for  the  electrodes  and  the  electrolyte 
involve  charge  conservation  in  the  solid  (Eq.  (7)  in  S-domain  in 
Fig.  5)  and  the  electrolyte  (Eq.  (8)  in  E-domain  in  Fig.  5)  and  mass 
transport  law  in  the  solid  (Eq.  (9)  in  S-domain  in  Fig.  5)  and  the 
electrolyte  (Eq.  (10)  in  E-domain  in  Fig.  5): 

V-(ffV0s)  =  0  (7) 


Fig.  3.  A  schematic  diagram  of  a  tilted  plane  (measured  plane)  and  horizontal  plane 
(correction  target  plane).  The  angle  between  the  two  planes  is  6  and  those  two  planes 
share  a  cross  line.  The  tilted  plane  is  rotated  about  the  cross  line  by  angle  6  to  correct 
the  image. 
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Fig.  5.  A  schematic  of  a  half-cell  and  model  representation.  Va  Geometry  was  employed 
using  symmetric  boundary  conditions.  ‘S’  and  ‘E’  represent  the  solid  and  the  electro¬ 
lyte,  respectively.  T  stands  for  interface  between  the  two  phases. 


nr  =  -Jp 


11,12 


-nr  =  -j\ 


I  SI  ,S2 


Kapp 


S3 


(12) 


For  the  charge  balance  in  the  electrolyte  expressed  by  Eq.  (8): 


nr  =  0 


I  El  —  E5 


nr 


=  Ji 


I II  ,12  ,E6 


(13) 


For  the  material  balance  in  the  solid  expressed  by  Eq.  (9): 


-n-r  =  0 


IS1 -S3 


,  -n-r  =  jY/F 


11,12 


(14) 


For  the  mass  transport  law  in  the  electrolyte  expressed  in  Eq. 

(10): 


-n-r  =  0 


IE1-E5 


,  -n-T  =  j]1  /F 


II  ,I2,E6 


(15) 


where  n  and  T  represent  the  normal  vector  and  flux  vector, 
respectively.  For  the  solid  mechanics  calculation,  roller  conditions 
were  applied  to  the  symmetric  plane  (SI  and  S2),  and  the  bottom 
surface  (S3)  was  fixed.  The  outer  surfaces  were  given  no  constraint. 
Uniform  initial  conditions  are  assumed,  i.e.,  ce  =  c°,  cs  =  c°, 
0 s  =  0s  *  0e  =  0e-  The  list  of  the  parameters  used  in  the  simulation 
is  provided  in  Table  1.  Solution  of  the  system  of  differential  equa¬ 
tions  was  performed  by  the  Finite  Element  Method  using  COMSOL 
software. 


4.  Results  and  discussion 


^  +  v(— DeVCe)  =  0  (10) 

where  a  is  the  conductivity  of  the  solid  particle,  0S  is  potential  in  the 
solid  particle,  k  is  ionic  conductivity  in  the  electrolyte,  0e  is  the 
potential  in  the  electrolyte, /is  the  mean  molar  activity  coefficient, 

is  the  cation  transference  number,  Ds  is  the  diffusion  coefficient 
in  the  solid,  Q  is  the  partial  molar  volume,  ah  is  hydrostatic  stress, 
De  is  the  diffusion  coefficient  in  the  electrolyte,  and  i  is  the  current 
density  in  the  electrolytic  solution. 

In  Eq.  (9),  the  first  term  in  the  flux  term  accounts  for  the  effect  of 
the  concentration  gradient,  where  lithium  ions  diffuse  to  lower 
concentration  regions  from  higher  concentration  regions.  The 
second  term  accounts  for  the  effect  of  the  stress  gradient.  This  effect 
originates  from  the  fact  that  an  intercalated  lithium  ion  displaces 
the  host  atoms  and  sets  up  a  strain  field  in  the  host  structure.  Since 
this  strain  field  affects  the  energy  of  a  second  intercalated  ion,  this 
leads  to  an  elastic  interaction  between  the  two  ions. 

The  stress-strain  relation,  including  the  effect  of  the 
intercalation-induced  strain,  is  given  by: 

1  rQ 

Eij  —  g  [(1  +  v)°ij  —  v(Tkk^ij]  +  <5y  (11) 

where  q/,  a#  are  strain  and  stress  components,  kk  represents 
summation  notation,  E  is  Young’s  modulus,  v  is  Poisson’s  ratio,  c  = 
c  -  c0  is  the  lithium  ion  concentration  change  from  the  stress-free 
state,  and  8y  is  the  Kronecker  delta.  The  last  term  accounts  for  the 
intercalation  strain.  Note  that  Eqs.  (9)  and  (11)  are  coupled  through 
the  concentration  c  and  the  hydrostatic  stress  ah. 

3.2.  Initial  boundary  conditions  and  material  properties 

The  equations  above  were  given  the  following  initial  boundary 
conditions: 

For  the  electronic  current  balance  in  the  solid  particle  expressed 
by  Eq.  (7): 


The  annealed  thin  films  of  LiCo02  were  analyzed  by  X-ray 
diffraction  using  a  PANalytical  X’Pert  Pro  system  with  a  molyb¬ 
denum  source  (A  =  0.709319  A).  Rietveld  refinement  was  done 
using  X-pert  High  Score  Plus  software.  The  calculated  lattice 
constants  were  a  =  b  =  2.8143  A  and  c  =  14.030  A,  which  were  in 
accordance  with  previously  reported  results  for  LiCo02  single 
crystals  [42]  and  synthesized  powders  [43].  Fig.  6  shows  the 
diffraction  pattern  of  the  thin  film  in  the  2-theta  range  containing 
the  three  principal  diffraction  lines  (003),  (101 ),  and  (104).  The  film 
showed  (104)  preferred  orientation  (46%).  With  the  (104)  plane 
parallel  to  the  substrate,  the  c-axis  makes  a  350°  angle  with  respect 
to  the  substrate  plane,  which  provides  access  to  the  lithium  layers 
for  faster  Li  ion  diffusion  as  compared  to  the  case  where  the  c-axis  is 
perpendicular  to  the  substrate. 

Fig.  7  shows  the  voltage  profile  of  a  cell  in  a  Swagelok  assembly 
with  a  constant  current  charge/discharge.  The  voltage  profile 
showed  a  plateau  at  around  3.9  V,  which  corresponds  to  the 
Li0.94CoO2  semiconductor  to  Li0.75CoO2  metal  transition.  This  data 
was  used  as  a  reference  for  the  in  situ  AFM  measurements.  One 
example  of  the  AFM  image  is  shown  at  the  top  left  in  Fig.  1,  where 
the  scanned  regions  are  at  the  center  of  the  figure,  and  the  two 
outer  flat  planes  represent  unscanned  regions.  These  AFM  images 
were  used  to  determine  the  height  changes  of  the  micro-pillars. 
Fig.  8  shows  the  average  heights  of  the  center  areas 
(340  nm  x  340  nm)  of  pillars  no.  3  (r0  =  3.8  pm)  and  no.  2 
(r0  =  2.05  pm)  shown  in  Fig.  1  while  the  applied  potential  was 
swept  at  a  rate  of  0.39  mV  s-1.  The  solid  blue  line  represents 
simulation  results  at  the  center  point  of  pillar  no.  3,  and  the  dotted 
red  line  represents  simulation  results  at  the  center  point  of  pillar 
no.  2.  The  blue  circle  bar  marked  ‘x’  indicates  the  experimental 
measurements  at  the  center  point  of  pillar  no.  3,  and  the  unmarked 
red  circle  bar  indicates  the  experimental  measurements  at  the 
center  point  of  pillar  no.  2.  As  seen  in  Fig.  1,  the  surfaces  of  the 
machined  pillars  show  roughness  due  to  multigrain  boundaries.  As 
Li-ions  are  removed  from  the  layered  structure,  the  volume 
expansion  will  vary  depending  on  the  grain  orientation.  The 
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Table  1 

Values  of  parameters  used  in  the  simulations. 


Parameters 

LiCo02 

Ref. 

LiCo02  particle  properties 

D  (mol  m  3) 

1  x  10  13 

[38] 

a  (S  m-1) 

10 

[38] 

E  (GPa) 

70 

[38] 

p  (kg  ITT3) 

5030 

[38] 

Q  (m3  mol ^) 

7.35  x  10  7 

[42,45-49] 

V 

0.2 

[38] 

Interface  interaction  properties 

Cathode 

«a, 

0.5,  0.5 

[50] 

k  (A  m-2- 

1.27  x  10"6 

[38] 

(mol  ITT3)"1-5) 
i'o  (A  ITT2) 

85.0  (anode  only) 

[51] 

Electrolyte  properties 

D2  (m2  s  3) 

7.5  x  10"11 

[50] 

k  (S  m-1) 

2.602  x  10“2  +  5.877  x  10“4ce  -  9.946x 
10_7c2  +  1.435  x  10-10c3  -  1.387  x  10“14c4 

[50] 

Constants 

R  (J-(molK)-1) 

8.314 

T(K) 

273.15  +  25 

F  (C  mol"1) 

96,487 

Initial  and  boundary  values 

c°(mol  m-3) 

1000 

d?(mol  m-3) 

0.98Cmax 

$(V) 

$(V) 

3.9  V 

0 

maximum  difference  can  be  as  much  as  50  nm  when  the  c-axis 
expands  by  2.4%  at  the  same  time  that  the  a-axis  shrinks  by  0.34%; 
if  the  initial  roughness  is  taken  into  consideration,  this  variance  can 
be  increased  further.  In  Fig.  8,  and  in  subsequent  figures  showing 
the  experimentally  obtained  height  change,  the  vertical  bars  serve 
to  reflect  the  surface  roughness.  The  half-range  of  these  bars 
represents  one  standard  deviation  from  the  mean  value  based  on 
the  measurements  of  height  change  over  the  scanned  area. 

In  order  to  observe  the  height  change  upon  delithiation  (charge) 
as  a  function  of  voltage  in  the  in  situ  AFM  cell,  the  average  value  is 
compared  to  the  simulation  result.  The  initial  voltage  in  the  simu¬ 
lation  was  set  to  3.9  V.  In  the  fluid  cell,  polarization  may  increase 
due  to  the  ohmic  resistance  across  the  electrolyte  where  the 
distance  between  the  reference  and  working  electrode  is  longer 


*  -  ai2o3 


8  10  12  14  16  18  20 

20  (deg) 


Fig.  6.  XRD  pattern  of  the  LiCo02  thin  film.  A1203  peaks  (*)  and  Au  peaks  (#)  originate 
from  the  substrate  and  current  collector  supporting  the  LiCo02  thin  film. 


Fig.  7.  Voltage  profiles  of  Li/LiCo02  from  a  thin  film  cell.  The  cell  was  cycled  between 
3.0  and  4.5  V  at  a  current  of  3  pA. 

than  in  the  Swagelok  case.  Also,  the  absence  of  an  external  pressure 
to  maintain  electrode-current  collector  contact  in  the  fluid  cell 
may  also  increase  the  contact  resistance  of  the  sample.  Thus,  under 
the  assumption  that  the  starting  points  of  the  simulation  are  the 
same  with  each  measurement,  the  relative  height  changes  from  the 
starting  point  are  compared.  While  scanning,  the  applied  voltage 
continued  to  increase,  so  the  middle  value  of  the  voltage  during 
each  scan  is  represented  in  the  resultant  figures  and  compared  with 
the  simulation  results. 

As  can  be  seen  in  Fig.  8,  in  the  case  of  pillar  no.  3,  both  the  height 
changes  calculated  in  the  simulation  and  those  found  by  experi¬ 
mental  measurement  show  similar  trends.  The  measured  height 
changes  were  found  to  gradually  increase  to  about  26  nm,  while  the 
simulation  calculated  a  change  of  22  nm  at  the  end  of  charge 
(4.15  V).  This  volume  expansion  is  caused  by  the  removal  of  the  Li- 
ions,  which  causes  adjacent  C0O2  layers  to  undergo  increased 
electrostatic  repulsions  as  negatively  charged  oxygen-oxygen 


3.6  3.7  3.8  3.9  4.0  4.1  4.2  4.3 

Potential  /  V 


Fig.  8.  The  height  evolution  of  pillar  no.  3  (r0  =  3.8  pm)  and  no.  2  (r0  =  2.05  pm)  when  the 
cell  was  swept  with  0.39  mV  s_1.  The  mean  value  of  the  height  over  340  nm  x  340  nm  as 
measured  experimentally  is  represented  by  markers,  and  the  vertical  bar  represents  the 
surface  roughness  as  one  standard  deviation  from  the  mean.  The  solid  line  shows  the 
simulation  result.  The  solid  blue  line  represents  the  simulation  results  on  pillar  no.  3, 
while  the  dotted  red  line  represents  the  simulation  results  on  pillar  no.  2.  The  blue  circle 
bar  marked  ‘x’  indicates  the  measurement  results  on  pillar  no.  3,  and  the  unmarked  red 
circle  bar  indicates  the  measurement  results  on  pillar  no.  2.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 
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interactions  increase.  It  is  worth  noting  that  the  abrupt  increases  of 
the  OCV  in  the  low  voltage  regions  result  in  slow  increases  in  the 
height  of  the  sample  because  of  low  current  flow.  However,  as  it 
passes  the  plateau  region  of  the  OCV,  the  flux  of  Li  ions  becomes 
increased,  even  at  constant  voltage  sweep  rates,  and  as  a  result,  the 
height  increases  quickly.  The  volume  increases  by  about  1.3%,  which 
is  lower  than  the  2.4%  expansion  in  Chex  from  LiCo02  to  in 
Lio.4gCo02.  This  is  because  the  grain  size  in  the  sample  is  50—70  nm, 
which  is  very  small  compared  to  the  pillar  size.  Thus,  the  volume 
expansion  can  be  considered  isotropic,  as  we  assumed  in  the 
simulations. 

Fig.  8  also  shows  results  obtained  from  the  small  pillar  (pillar  no. 
2  in  Fig.  1,  ro  =  2.05  pm).  Similar  to  the  case  of  the  larger  pillar  (pillar 
no.  3  in  Fig.  1,  ro  =  3.8  pm),  the  simulation  results  follow  the 
measurement  very  closely.  However,  here  the  maximum  height 
change  is  smaller  than  in  the  case  of  the  larger  pillar.  The  height 
increase  comes  from  volume  expansion  resulting  from  electrostatic 
repulsion.  It  implies  that  a  larger  pillar  will  exhibit  a  larger  height 
change.  The  effect  of  pillar  size  on  the  deintercalation-induced 
vertical  displacement  is  demonstrated  in  Fig.  9,  where  the  height 
change  as  a  function  of  potential  is  compared  for  the  two  pillars 
with  different  radii.  It  can  be  seen  that  during  the  initial  stages,  at 
the  same  current,  the  height  change  in  the  smaller  pillar  is  larger 
because  Li-ion  extraction  occurs  not  only  from  the  surface  but  from 
the  inner  portion  of  the  pillar  as  well;  this  is  due  to  the  smaller  size 
of  the  pillar.  In  contrast,  in  the  larger  pillar,  Li-ions  are  extracted 
from  the  surface,  resulting  in  lower  expansion.  As  the  cell  is  charged 
further,  however,  the  larger  pillar  also  experiences  similar  dein¬ 
tercalation  at  its  interior,  resulting  in  greater  height  change. 

The  sweep  rate  is  another  important  factor  in  determining  strain 
and  stress  due  to  intercalation  because  it  determines  the  Li-ion 
concentration  distribution  inside  the  particle.  The  current  density 
at  the  surface  of  a  particle  inside  a  cell  may  vary  depending  on  the 
C-rate  and  also  on  the  location  of  the  particle  in  the  electrode;  the 
current  density  is  higher  for  particles  near  the  separator  (electrode 
surface).  To  elucidate  the  influence  of  the  sweep  rate  on  the  volume 
strain  level,  the  sweep  rate  was  increased  by  a  factor  of  two. 
Contamination  from  the  side  reactions  and  gas  generation  resulted 
in  a  low  quality  AFM  image.  Therefore,  the  discussion  of  the  de- 
lithiation  rate  effect  is  based  on  the  results  of  simulations 
(Fig.  10).  The  measured  global  height  changes  at  the  higher  sweep 
rate  are  smaller  than  those  measured  at  the  lower  sweep  rate.  As 
the  sweep  rate  increases,  the  total  amount  of  Li-ion  extracted  from 
the  particle  decreases  so  that  the  volume  expansion  also  decreases. 
However,  a  higher  current  density  may  induce  a  higher  gradient  in 


Fig.  9.  The  numerical  results  of  the  micro-pillar  size  effect  on  height  evolution.  The 
cells  were  swept  with  0.39  mV  s-1  rate. 


Fig.  10.  The  numerical  results  of  the  sweep  rate  effect  on  height  evolution.  The  cell  was 
swept  with  0.39  mV  s-1  and  0.78  mV  s-1  rates. 


the  lithium  concentration,  as  can  be  seen  in  the  results  of 
a  numerical  simulation  shown  in  Fig.  11.  Here,  pillar  no.  3  was  swept 
with  different  sweep  rates  (a)  0.39  mV  s-1  (b)  0.78  mV  s-1.  This 
high  gradient  can  lead  to  a  higher  stress  level.  Also,  a  high  discharge 
rate  may  cause  an  unintended  over-discharge,  which  may  lead  to 
degradation  of  the  cell;  if  the  current  is  sufficiently  high,  the  state  of 
charge  (SOC)  of  the  material  could  exceed  x  =  1  at  the  sharp  surface 
of  particles  in  the  electrode. 

Numerical  modeling  is  an  advantageous  approach  for  evaluating 
the  effect  of  the  different  operating  conditions  on  strain  and  stress 
levels,  whereas  it  is  not  easy  to  measure  stress  level  via  an  AFM 
experiment.  The  current  modeling  results  are  consistent  with 
measurement  results  in  principle.  This  implies  that  numerical 
results  can  provide  further  information  than  was  measured. 
Fig.  12(a)  shows  the  calculated  displacement  in  the  z  direction  and 
Fig.  12(b)  shows  total  displacement  in  pillar  no.  3  when  it  was 
swept  with  0.39  mV  s_1  rate.  The  deformation  is  exaggerated  100 
times  to  show  the  shape  change.  Due  to  the  unconstrained 
boundary  condition  and  lower  Li-ion  concentration  around  the 
outer  surface,  the  specimen  expands  further  at  the  outer  surface. 
Consequently,  the  pillar  shape  is  convex.  This  convex  shape  is  also 
verified  from  the  measurements  shown  in  Fig.  13  for  pillar  no.  3. 
The  average  height  is  decreased  as  the  radial  position  is  increased. 

Yield  stress  and  von  Mises  stress  are  not  proper  failure  criteria 
for  brittle  materials.  Also,  the  compressive  strengths  for  ceramics 
are  much  larger  than  the  ultimate  tensile  strength.  The  maximum 
tensile  stress  is  appropriate  information  in  predicting  fracture 


Fig.  11.  The  lithium  concentration  profile  in  pillar  no.  3  when  the  potential  is  4.15  V 
with  different  sweep  rates:  (a)  0.39  mV  s-1  (b)  0.78  mV  s-1. 
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a  b  (nm) 


Fig.  12.  Exaggerated  deformation  of  pillar  no.  3  where  the  color  map  represents  (a) 
z-displacement  (b)  total  displacement. 


failure  of  brittle  materials.  The  maximum  principal  stress  occurs  at 
the  center  of  the  micro-pillar  due  to  expansion  of  the  specimen. 
Fig.  14  shows  the  computed  maximum  principal  stress  in  the  micro¬ 
pillar  as  a  function  of  potential  for  different  sweep  rates  and 
different  applied  boundary  conditions.  As  the  cell  is  charged,  the 
maximum  stress  increases  to  around  217  MPa  when  the  voltage  is 
4.3  V  and  the  boundary  of  the  bottom  of  the  pillar  is  fixed.  However, 
if  this  boundary  condition  is  changed  to  a  roller  condition 
(unconstrained  in  the  x—y  plane,  but  with  zero  displacement 
prescribed  in  the  z-direction),  the  stress  level  is  decreased  to 
35  MPa.  This  is  attributed  to  a  relaxation  of  the  strain  by  allowing 
free  movement  in  the  horizontal  plane.  For  the  same  reason,  the 
vertical  displacement  is  larger  with  a  fixed  boundary  condition. 
Similarly,  the  current  density  effects  are  different  depending  on  the 
boundary  condition.  For  a  roller  boundary  condition,  higher  current 
density  indicates  a  higher  stress  level.  However,  for  the  fixed 
boundary  condition,  the  stress  level  is  not  changed  significantly, 
and  a  lower  sweep  rate  results  in  a  slightly  higher  level.  This  is 
attributed  to  a  relaxation  in  the  roller  boundary  condition,  so  the 
concentration  gradient  is  more  critical  than  the  actual  Li- 
concentration.  This  can  also  explain  why  the  maximum  stress 
level  occurs  during  the  charge  process  where  the  concentration 
gradient  is  maximized.  However,  in  the  fixed  boundary  condition, 
there  is  no  relaxation  phenomenon,  which  means  that  as  Li-ions  are 
extracted  further,  there  is  more  volume  expansion.  This  causes  the 
monotonous  increase  in  stress  level  in  the  time  profile. 

The  deviation  in  the  measured  height  change  from  the  simula¬ 
tion  may  result  from  many  different  factors  in  addition  to  the 


Fig.  13.  Comparison  between  experimentally  measured  results  and  those  predicted  by 
simulation  based  on  cross-sectional  profiles  of  the  micro-machined  pillar  no.  3  after 
de-lithiation. 


Fig.  14.  Mechanical  boundary  condition  effects  with  different  sweep  rate  on  the 
maximum  principal  stresses  in  pillar  no.  3. 


intrinsic  measurement  error  associated  with  the  instrumentation 
and  technique.  LiCo02  is  known  to  dissolve  into  the  electrolyte, 
which  may  impact  the  measured  heights  of  the  electrode;  weights 
of  dissolved  metal  ions  measured  by  Atomic  Absorption  Spectros¬ 
copy  (AAS)  have  been  reported  for  LiCoCb  as  0.8%  of  the  sample 
weight  [44].  Also,  the  adhesion  of  the  active  material  to  the 
substrate  may  become  worse.  This  phenomenon  was  observed  in 
a  different  sample.  After  continuous  scanning  in  contact  mode  with 
a  Si3N4  tip,  fracture  occurred  at  the  point  of  attachment  of  the 
micro-machined  pillar  to  the  current  collector.  Adhesion  between 
the  gold  substrate  and  the  LiCoCb  film  was  not  sufficient  to  endure 
continuous  force  exertion.  Also,  as  there  was  gas  generation  during 
cycling,  the  resultant  bubbles  hindered  scanning  and  prevented 
obtaining  good  images.  In  addition,  the  AFM  probe  shielding  effect 
may  have  changed  the  Li-ion  content  over  the  pillar  sample. 
Another  characteristic  feature  of  LiCo02  that  might  have  affected 
the  results  is  phase  transition;  X-ray  diffraction  examinations  have 
indicated  that  a  two  phase  region  can  be  generated  between 
x  =  0.75  and  x  =  0.975  for  a  LixCo02  electrode  [14]. 

5.  Conclusions 

This  work  used  a  scanning  probe  technique  to  understand 
dynamic  deformation  processes  during  Li  deintercalation  in  a  Lix_ 
C0O2  cathode  material.  A  unique  experimental  approach  with 
a  pre-defined  micro-machined  specimen  has  allowed  for  in  situ 
characterization  of  morphology  of  the  specimen.  The  dimensional 
changes  of  FIB  micro-machined  samples  were  monitored  in  situ, 
and  the  results  showed  expansion  in  volume  due  to  electrostatic 
repulsion.  The  volume  expansions  are  generally  in  good  agreement 
with  previous  XRD  measurements.  The  averaged  height  increase 
calculated  from  XRD  measurements  of  volume  expansion  is  about 
1.34%  when  x  is  changed  from  1  to  0.35  in  LixCo02,  while  the 
measurement  from  this  in  situ  AFM  experiment  is  about  1.28-1.3%. 
The  results  of  the  modeling  follow  the  experimental  results  rather 
well.  The  changes  in  volume  expansion  due  to  the  geometry  and 
kinetic  rate  changes  showed  the  same  trend  with  our  simulation 
results— a  steep  expansion  in  plateau  OCV  regions  and  a  larger 
expansion  for  the  larger  pillar.  The  computed  profile  of  the 
deformed  by  de-lithiation  micro-pillar  matches  well  with  the 
experimentally  measured  shape  changes  in  the  specimen.  Further, 
the  simulation  results,  validated  from  measurements,  provided  the 
stress  state,  which  can  be  used  in  assessing  the  structural  stability 
of  the  material  but  is  not  easily  measured. 
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